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Abstract. Reliable single photon sources constitute the basis of schemes for quantum 
communication and measurement based quantum computing. Solid state single 
photon sources based on quantum dots are convenient and versatile but the electronic 
transitions that generate the photons are subject to interactions with lattice vibrations. 
Using a microscopic model of clcctron-phonon interactions and a quantum master 
equation, we here examine phonon-induced decoherence and assess its impact on 
the rate of production, and indistinguishability, of single photons emitted from an 
optically driven quantum dot system. We find that, above a certain threshold of 
desired indistinguishability, it is possible to mitigate the deleterious effects of phonons 
by exploiting a three-level Raman process for photon production. 
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1. Introduction 

Single photon sources are an essential component of many quantum information 
processing (QIP) protocols, from quantum key distribution (QKD) protocols [U E] to 
linear optical quantum computing (LOQC) schemes [3J IH E]. Optical schemes using 
path erasure, a two-photon interference effect, can be used to generate long-range 
entanglement between physically separated systems [6j IH |9]. Such procedures can 
be repeated on many different pairs of systems and so create a distributed cluster 
state pm E], which is the key resource required for implementation of measurement 
based quantum computing. 

In order to be useful in these applications, a photon source must be of a high quality 
in two respects: it must reliably produce a single photon on demand, and the photons 
produced must be indistinguishable from one another. 

To be perfectly indistinguishable, photons must have the same pulse width, 
band width, polarization, arrival time at the detector, and carrier frequency. 
Indistinguishability is vital if photons are to exhibit high quality quantum interference, 
and its consideration is therefore critical when designing photon sources for LOQC 
and path erasure entanglement generation. If photons can be distinguished even in 
principle, this can lead directly to imperfect LOQC gates, or to a lessening of the degree 
of entanglement generated in distributed cluster states. In QKD, indistinguishability is 
less important as interference effects are not required, but it is of paramount importance 
that no more than one photon is emitted on demand; multiple photon emission leads to 
security loopholes [5]. 

As a measure of indistinguishability we exploit the Hong-Ou-Mandel (HOM) effect 
[TT] . which relies on the bunching behaviour of identical photons when they are incident 
on a beam splitter with the same temporal profile. If the photons are indistinguishable, 
they will always emerge in and be detected in the same output arm (Fig. [I]). The 
number of same-arm detection events are usually plotted as a function of arrival time 
of the two photons - and hence a 'dip' at a time difference of zero is an indicator of 
indistinguishability. The degree of distinguishability is typically given by the visibility 
of the dip, f hom, which is the proportion of same-arm detections, over many runs of the 
experiment. 

As mentioned above, besides indistinguishability an essential characteristic of a 
good photon source is that it will consistently produce a single photon, but never more 
than one, on demand [12]. A laser source can be attenuated so that it gives zero 
or single photons most of the time. However, photons from classical lasers sources 
obey Poissonian statistics so that in order to keep the two-photon event rate low the 
probability of a single photon may become unworkably small for many applications. 
In addition to this, Poissonian sources are unsuitable for two-photon interference 
experiments since the rate of two photon production from a single source is similar to 
that for a single photon from each source [5]. It is worth noting that perfect efficiency 
should never be a requirement in any realistic optical QIP scheme, as these schemes 
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Figure 1. Hong-Ou-Mandel effect: a pair of photons from two sources S\ and 
S2 incident on a beam splitter exhibit perfect bunching behaviour if they are 
completely indistinguishable. In this case the detectors D + and D_ will never click 
simulataneously. 

must always be tolerant to photon loss within other parts of the apparatus. However, 
a good photon source must be reasonably efficient to be useful, and we should not be 
forced to trade efficiency for other desirable characteristics. 

The use of low-dimensional quantum systems as single photon sources avoids the 
efficiency problems of Poissonian sources. Successful experimental implementations 
have been realised in a number of different systems including atom-cavity schemes 
[USUI EE HE!, quantum dots [HI EE] and diamond colour centres [HE US [21]. While the 
majority of the work has been focussed on the efficient production of a single photon, 
recent experiments in NV centres [20j [21] and quantum dots [22j [23] have demonstrated 
two-photon interference effects from different sources, albeit sacrificing efficiency by 
filitering out undesired frequencies. It has been suggested that cavities could be used 
in these systems, to enhance the emission into the target mode, reducing the need for 
filtering [23] . 

In order to improve the characteristics of a photon source, it is not sufficient to 
simply consider the material parameters of the system being used: one should also 
consider the approach used to control the system. Perhaps the simplest strategy is to 
excite the system first optically, either coherently or incoherently, and wait for the system 
to relax into its ground state, emitting a photon in the process; we will henceforth refer 
to this as the 'pulse-relax' technique. This approach makes minimal resource demands 
on the system and, due to its simplicity, is the technique proposed in some remote 
entanglement generation schemes [61 [8]. The pulse- relax approach is problematic in 
systems where the excited state is sensitive to decoherence, which will degrade the 
photon's indistinguishability [251 ESI El [27]. These effects can be reduced, for example 
by exploiting the Purcell effect to enhance the emission rate into the desired photon 
mode [2H1 [29] or by using temporal post-selection of emitted photons [27]. However 
with experimental limits on cavity couplings both of these inevitably lead to lower 
efficiency as the proportion of emissions that are utilized falls [30] . 

A fundamentally different approach is to use more elaborate QED schemes to release 
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a photon from the system in a controlled manner. In particular single photon sources 
using a Raman approach have been analysed [25] [26] and experimentally realised [13] . 
The approach places more demands on the system, requiring a three level system with 
a A-system configuration. Each arm is coupled to either a classical or a quantum light 




Figure 2. Effective level structure of a driven A— system inside an optical cavity. 
Here, fi is the amplitude of an external (laser) driving field, h the optical dipole-cavity- 
coupling, and v the shared detuning of the laser and the cavity transition frequency 
from the excited state. The system can decay into state |<7i) and in doing so emit a 
photon into a well-defined external mode. 

field; in Fig [2] we show the situation when one arm is driven classically, with a coupling 
Q, and the other arm is coupled to a cavity mode with strength g. By detuning both 
arms of the A-system by the same amount v, population is transferred from one arm to 
the other, whilst suppressing population in the top state. Provided that the coupling 
strengths are small in comparison to the detuning, Q <C v and g <C u, we induce an 
effective coupling between the two low-lying states, which causes oscillations with Rabi 
frequency gQ/Au. The population in the excited state remains small at all times, and 
so any decoherence that arises due to environmental coupling to this excited state may 
be reduced using this strategy. In particular, a common realisation of a single photon 
source is a quantum dot, in which the excited state typically has a different charge 
configuration to the ground state, and so couples to acoustic phonons (31] that can act 
as a noise source. 

In this paper we provide a detailed and realistic analysis of the effects of the lattice 
vibrations for single photon emitters based on (self-assembled) semiconductor quantum 
dots. By comparing a standard pulse-relax approach with the aforementioned Raman 
technique, we find that the latter can can offer vast improvements in terms of both 
photon indistinguishability and source efficiency. Our results extend and complement 
a previous study [26] which considered generic pure dephasing noise. However, here we 
also show that the precise choice of control parameters is important if phonon-induced 
decoherence is to be successfully suppressed. 
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2. Our model 

For convenience of notation we consider the A-system detailed in Fig. [2] for both the 
pulse-relax and the Raman approach. In both cases one arm of the system is coupled 
to the cavity with detuning v. In the Raman scheme, the other arm is driven with 
strength Q by a laser with a matching detuning u, and the system starts in state \go). 
By contrast, we model the pulse-relax approach by setting Q = v = and starting in 
state |e), ignoring the details of the excitation process. The fact that we are neglecting 
the excitation step in this simplified picture will slightly favour the pulse-relax approach, 
but is largely justified on the assumption that the initial excitation process takes place 
quickly compared with other system dynamics. Our framework thus allows us to consider 
the pulse-relax approach as a special case of the master equation we will now derive for 
the Raman approach. 

We split the Hamiltonian into contributions from the emitter system and cavity 
{sys), the driving laser (dr), the unperturbed phonon bath (ph) and a system-phonon 
interaction term (int): 

H = H sys + H dr + H ph + H int ■ (1) 

In the following, we shall describe each of these parts separately. 

Using the standard Jaynes-Cummings Hamiltonian, the A-system coupled to a 
cavity mode with strength h and detuning v from resonance can be described by the 
following Hamiltonian: 

H sys = (u + v) \e)(e\ + ua)a + ^ (|pi)(e|a t + h.c.) , (2) 

where h.c. denotes the Hermitian conjugate, u is the cavity mode frequency. 

We can restrict ourselves to one or zero cavity photons since the Jaynes-Cummings 
model preserves excitation number and, due to the structure of the system, once a 
photon has escaped re-excitation is impossible. This allows us to replace \gi) a) with a 
new combined atom-photon state \g a ) and a) a with |<7 a )(<7 a |. 

The term H dr describes a laser driving the transition \qq) -h- |e) with coupling 
strength Q, and whose frequency is detuned from resonance by the same v parameter, 
yielding 

H dr = fJcosfwt) (|e}(# | + h.c.) . (3) 

After performing the rotating wave approximation, assuming that u 3> z/, h, Q, we are 
left with 

H sys + H dr = u\e)(e\ + l - (h\g a )(e\ + Q\e)(go\) + h.c. (4) 

for the total system Hamiltonian. 

We model the phonons as a bath of harmonic oscillators 

H ph = ^2^qblb q , (5) 
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which couple to the exciton state |e) through the deformation coupling with coupling 
constants f q = D\q\ in the usual way [32J: 

H int = \e)(e\Y,f q {b\ + b q ). (6) 
q 

The effect of the phonon bath on a driven quantum dot can be described using 
a Lindblad master equation [33], where the Lindblad operators induce phonon-assisted 
transitions between the dressed system eigenstates [31] . Being analogous to Refs (33J [31] , 
it suffices to outline the derivation of the phonon master equation only briefly in 
the following. To derive the master equation we must first diagonalise the system 
Hamiltonian, H sys + Hd r - We find that the eigenvalues are 

A = 0, 

v ± y/u 2 + n 2 + h 2 
A± 2 , (7) 

with corresponding eigenvectors 

IV'o) = «o (ft \go) - fi \g a )) , 

W±) = n± (ft \g ) + h \g a ) + 2A± |e)) . (8) 

uq and n± are appropriate normalisation factors. Making the Born and Markov 
approximations leads to the master equation [33] in standard Lindblad form: 

p = t[p,H}+D ph (p) (9) 

with the phonon dissipator given by 

D ph (p) = J(A) [(iV(A) + 1)D [Pa] p + N(A)D[P A ]p\ (fO) 

where D[L]p = LpV — l/2(L^Lp+pL jf L), A = A+-A_, and P A = -|^_)(^+|. Note that 
the phonons only induce transitions between the two optically bright system eigenstates 
and do not couple to the dark l^o)- m the above equation N(A) is the bosonic mode 
occupation number: 

N M = ^rr (") 

P — {kBT)~ l and we shall be henceforth consider all systems at room temperature, 
T = 298 K. The spectral density function J(u) represents the electron-phonon coupling 
weighted by the density of phonon modes [35]. We expect this to be dominated by 
deformation potential coupling, and in this case we obtain [33] : 

J(u) = «u 3 e"^) 2 . (12) 

We take a = 0.0027 ps _1 and co c = 2.2 ps _1 , values that agree well with experiments on 
self-assembled quantum dots J3TJ [3B] . 

We absorb the rates in Equation ( 10 ) into the Lindblad operators to obtain following 
decoherence operators: 



U + = y/J(N+l)W-)(iP+\ (13) 
U„ = yfjN\il) + ){^\ (14) 
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where we have taken J = J (A) and N = iV(A). 

As a measure of the degree of indistinguishability of the photons produced in the 
emission process, we consider the HOM visibility, which is the normalised probability 
of same arm detections obtained over many runs of the experiment, 

fHOM = : , (15) 

Psame + PdiS 

where p samc = p(D+\D+) + p(D_\DJ) and p diS = p(D + \D-) +p(D-\D + ) with p(D x \D y ) 
being the probability of obtaining a click in detector D x conditional on the previous 
click having occurred in detector D y . 

In order to calculate whom we thus need to consider photons emitted from two copies 
of the system Si and S2, one for each input arm of the beam splitter interferometer. 
The joint state of the system then inhabits the space S = Si ® S2. We label the two 
beam-splitter input modes ai and a 2 , and the two detection modes c + and c_. The 
beam splitter performs the transformation 

K—=(ai + a 2 ) (16) 
V 2 

k—= (ai - a 2 ) (17) 
v2 

where k is the cavity leakage rate. We assume that the field in the mode outside the 
cavity is directly related to the field inside, neglecting the process of escape from the 
cavity. The effect on the system of a detection in the plus or minus output mode is then 
described by 

C ± = \gi){g a \®l±l®\gi){g a \. (18) 

We could now simulate many trajectories of this system and build up an estimate of 
whom by averaging these [37] • Instead we use a semi-quantum master equation technique 
described in more detail in Appendix A This allows us to find vhom in a single run of 



a master equation acting on a slightly larger Hilbert space. 




Figure 3. The system jump-space: On the first jump the system moves to P+ /P_ 
depending on the arm in which the photon is detected. After the second jump the 
system moves to state Ps/Pd depending on whether the second photon was detected 
in the same or different arm to the first. At any point the system can undesirably 
spontaneously emit into the environment, moving to the junk state Pe- 
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When calculating vhom we must consider events in which photons trigger either 
detector, and events in which photons are spontaneously emitted into the environment. 
We introduce a set of process-states to record and labels these event classes (Fig. [3]): 

S P = {P , P+, P^,P S , P D , Pe} ■ (19) 

The process-state starts as Pq and remains there until an event of interest occurs. P + /_ 
represents the process-state after a single photon has been detected in the -D+/- detectors 
respectively. After a second photon has been detected the process-state becomes Ps/d 
depending on whether the second photon was detected in the same or different detector 
as the first. If at any point a photon is emitted into the environment the process moves to 
state Pe- When calculating the indistinguishability we can then ignore any population 
in state Pe, but we must include it when considering the overall efficiency of the process. 

We must also identify the operators that cause the movement between the process- 
states. Before doing this we extend the state-space S of the system to include the 
process-states: 

S = Sx ® S 2 <8> S P . (20) 
The detection operators are then given by 



a 
c. 
a 
c. 



(0) 

f 

(+) 
(-) 

f 

(0) 



C+ <g> 
C+ <g> 
c+ <g> 

c_ ® 



\P+){P*\ 

\Ps)(P+\ 
\Pd)(P-\ 
\P-)(Po\ 
\Ps)(P-\ 
\Pd)(P+\ 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 



where for example ^ is the jump operator representing a second detection in the D + 
detector, when the first detection was in D_. The spontaneous emission operators are 
given similarly: 

(27) 
(28) 



E. 
E 
E. 
E 
E. 
<+) 



(0) 


= \g 2 )(e\ ®l5® 




(0) 
2 


= h® \g 2 )(e\ ® 


l^)(Po| 


(+) 


= \g 2 )(e\ 


l^)(P+ 


(+) 
2 


= I s ®\g2)(e\ ® 




(-) 


= \92)(e\®I s ® 




(-) 
2 


= I S <8> \g 2 )(e\ ® 





(29) 
(30) 
(31) 
(32) 



where here E\ represents a emission from Si acting after the first photon was detected 
in the D + detector. 

Finally we must modify our phonon decoherence operators. The decoherence 
process must happen independently on each subspace, as classically separated branches 
of the process are unable to interfere with one another. 



U 



(0) 



VJ(JV + l)|V_)(V'+l®Is ® \Po)(P \ 



(33) 
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with similar operators acting on the second system. We do not need decoherence 
operators acting on the Ps, Pd or Pe, since we are only concerned with populations in, 
and not coherences between, these states. Moreover, we only really need to keep track 
of the total population in each of these subspaces - a fact that we exploit to reduce the 
dimension of our problem for the numerical simulations. 

We form a Lindblad master equation using these 24 Lindblad operators: 

p = t[p,H] + Yl H { L iP L \ ~ l / 2 ( L l L iP + P L \Li)) ■ (39) 

i 

The ji are the rates for each process. As noted earlier for the U± \ this rate is 1, as the 
rates have been encorporated into the Lindblad operators. For the C± we need 7 = k, 
the cavity leakage rate, which we take to be 3 * h. For the we take 7 = 0.05ps -1 , 
assuming a ratiative livetime of 200 ps. 

The dimension of this extended space S is 4x4x6 = 96. In fact we can reduce this 
by eliminating some unneccessary states from the subspaces. By carefully considering 
the basis states accessible in each the subspace corresponding to each process, we can 
reduce the number of system states to 9 + 6 + 6 + 1 + 1 + 1 = 24. For our simulation 
we will need to calculate the density matrix for this system. As noted in |Appendix A 
no coherences can exist between the different process-state subspaces. This reduces the 
number of density matrix elements we need to track to 9 2 + 6 2 + 6 2 + 1 + 1 + 1 = 116. 

3. Results 

Fig. [4] shows the HOM visibility obtained from, and spectral density used in, simulations 
of the pulse-relax technique. At low coupling strengths the phonon spectral density is 
small, and so phonon decoherence is largely avoided giving high indistinguishability. 
However, in the extraction of the photon into the cavity mode is slow and thus 
spontaneous emission into environmental optical modes is a problem. To take account 
of this effect, we define the 'combined HOM visibility', which is the product of success 
probability and bare HOM visibility, which approaches zero as h — > 0. Taking large 
coupling strengths allows one to access the region to the higher frequency side of the 
hump in the phonon spectral density and so could avoid this problem, but is unrealistic 
given the cavity parameters currently obtainable experimentally. In the experimentally 
feasible region (h < 1) [3H] we must therefore trade indistinguishability for efficiency. 

In contrast, the Raman procedure (Fig. [5j) avoids this trade-off. For small detuning 
the visibility is low, but this is because we do not get a proper Raman ground state 
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cavity coupling strength, h (ps ') 



Figure 4. Pulse-relax technique: calculations are performed as a function of the cavity 
coupling strength. Upper panel: at realistic coupling strengths (h < 2) increased HOM 
indistinguishability necessarily entails a decrease in efficiency (success probability), 
with the product of the two (blue curve) approaching zero. Here the spontaneous 
emission rate is T — 0.05 ps _1 . For overcoming phonon-induced decoherence and 
achieving a high success probability, we must move into a region of unreasonably high 
h. The lower panel shows the phonon spectral density, Eq. [T2J evaluated at the cavity 
coupling strength h, giving a rate that is directly proportional to phonon-induced 
dephasing during the pulse-relax process. 




Figure 5. Raman technique, h = 0.5: by choosing a detuning to move beyond the 
region of high phonon spectral density we can achieve near-perfect indistinguishability. 
The efficiency in this region is high enough for a feasible photon source. 



transition unless h v. If this condition is not met, the system simply undergoes 
a non-optimal detuned pulse-relax transition. Once we reach a detuning of around 
12h = 6 the indistinguishability and efficiency both increase. Our choice of detuning 
size is limited only in that in must be small in comparison with the original energy gap 
between the ground and excited states. This leaves us free to use large detunings to 
push to frequencies above the region of high phonon spectral density. As the detuning 
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is increased the efficiency saturates below unity. In this region both the time taken for 
the process and the average lifetime scale with the detuning squared, the two effects 
cancelling one another. 




0.0010 




10 D 10"* 10" 3 
photon indistinguishability 



Figure 6. Rate of photon pair production: with no spontaneous emission (upper), 
and spontaneous emission with an excited lifetime of 200 ps (lower) 

In many applications where photons are used, the product of efficiency and 
indistinguishability may not be the most useful metric for characterising the 
performance of the source; often photon escape errors can be accounted for, and 
the indistinguishability of photons that are detected is the important figure of merit. 
We therefore also calculate the rate of production of pairs of photons of a given 
indistinguishability using each approach (Fig. [6]). For our production rate we take 

Tf = ^ (40) 

where e/ is the efficiency and tf is the time taken for 99% of the runs to have completed 
(possibly unsuccessfully), for parameters chosen to obtain a given indistinguishability, /. 
This figure is somewhat approximate as it takes no account of how successful runs (where 
two photons are emitted into the correct modes) and unsuccessful runs are distributed 
within the process run time, and no allowance is made for time taken to reset the system 
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in the event of a failure. The effect of the former is minor since in our model spontaneous 
emission can occur uniformly at any point of the process. We will revisit the effect of 
the latter shortly. 

Even in the absence of spontaneous emission (Fig. |6j upper panel), the Raman 
procedure is quicker than the pulse-relax process at generating photons of a sufficiently 
high level of indistinguishability. With the dephasing parameters chosen in our model 
this occurs for indistinguishability of greater than 99.99%. In our model, spontaneous 
emission is the only process degrading the efficiency - without it we have perfect 
efficiency and so neither of the potential shortcomings discussed in the previous 
paragraph apply. 

When spontaneous emission is added, we see a similar pattern but the 
indistinguishability threshold is lower (99.9% working with a spontaneous emission rate 
of 0.05 ps _1 ). This is an upper bound, as here the reset time becomes important. The 
efficiency of the Raman procedure remains fixed at about 80%, requiring on average 
1.25 runs per pair. In contrast, the efficiency of the pulse- relax procedure heads towards 
zero, meaning that many attempts will be needed to produce a pair. If the time taken 
to reset the system (to the excited state |e) in which we have assumed the pulse-relax 
system starts) is large, the Raman procedure will become advantageous at a far lower 
threshold. 

4. Conclusion 

To conclude, we have developed a realistic and microscopically justified model of the 
impact of phonons on solid state single photon sources. We used a modified, 'semi- 
quantum', master equation method for the efficient calculation of coincidence rates, 
without having to resort to a quantum Monte-Carlo simulation approach. In systems 
where phonon dephasing and spontaneous emission are the dominant loss channels, and 
with current limits on the obtainable cavity coupling h, we find that a Raman technique 
is preferable to the pulse-relax approach for producing highly indistinguishable photons 
that would be suitable for the most demanding applications of quantum information 
processing. 
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Appendix A. Semi-quantum Master Equations 

We are often interested in 'observable events' in quantum systems, such as the emission 
of a photon. When an event is observed the system undergoes a transition due to 
wave function collapse. During a period when no event is observed the system evolves 
according to a conditional master equation, reflecting the fact if events could be observed 
but are not, then this also informs our knowledge of the state. In order to answer 
questions about the probabilities and time distributions of events or chains of events, one 
approach is to use a quantum jump master equation to generate individual trajectories 
of the system. In each timestep we decide probabilistically whether an event should 
occur. If it does occur then the system collapses according to a quantum jump; if it 
does not occur then the system evolves conditioned on no jump occurring. Statistics 
about the quantities of interest are built up as many trajectories are created. 

Here we look at a different approach to calculating the quantities relating to events 
that occur in such systems. Instead of simulating multiple trajectories of the system, 
we efficiently increase the size of the statespace to record the information of interest. 
This allows us to calculate the desired system properties, and their time evolution, by 
solving a single master equation. 

Consider a system that can exist in a number of different states. Let a movement 
between these states constitute an event, and assume that each kind of event happens at 
a given rate. Such a system is heavily reminiscent of a classical continuous time Markov 
chain (CTMC), which can be represented as a graph with the states as nodes and the 



edges events weighted by the transition rates (Fig. Al). Given an initial state i in a 




Figure Al. A classical continuous-time Markov chain, with statespace W = 
{1, 2, 3, 4, 5}. The edge weight, r«, represents the transition rate from state i to state 
j. If Pi(t) is the population in state i at time t, the system is governed by the rate 
equations p\ = V,. u r ,,,>,. 



chain of size n, we can calculate the probability that at a later time t the chain is in 
state j, by solving the rate equations - a set of n ordinary differential equations. 
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Quantum systems differ from classical continuous-time Markov chains due to 
quantum superposition. We are not able to simply record the population in each 
quantum state as inter-state coherences are also important. Using a Markov chain 
to model the whole quantum system is not possible by definition - systems that can be 
modelled in this way do not exhibit quantum behaviour. 

In what follows it is helpful to explain carefully what we mean by 'state'. In 
quantum systems a state is usually a vector in the Hilbert space of the system. We 
shall call this a quantum-state. We can also refer to the 'state' of the overall process, 
considering for example a system that has emitted a photon to be in a different process- 
state to one which has not. A system changes process-state when an event is observed. 

Whereas thinking in terms of the CTMC is not useful when considering quantum- 
states, it is an effective way to think about process-states. As process-states are 
separated by an observed event, no coherences can exist between the two histories, 
making a CTMC approach feasible. Of course, process-states alone are not enough to 
model the whole system. Each process-state needs its own copy of the system attached 
to it. We can think of a Markov chain with a copy of our system at each node, where 
the transition rates are determined by the jump Lindblad operators corresponding to 
the events. Another way of thinking about this is that we extend our overall space with 
a set of process-states, to allow us to record events in the system. 

Formally, we take a set of process-states Sp, transitions between which correspond 
to our observable jump events described by jump operators Jq\ We extend the Hilbert 
space of our quantum system Sq by forming the tensor product: 

S = S Q ®S P . (A.l) 

The new Hamiltonian is given by 

H = H Q ®I. (A.2) 

If event Jq 1 causes a transition from system state a to b we say it is of type (a, b). Its 
action on the extended system S is described by 

J(0 = jW <g, |6>< a |. (A.3) 

For any other Lindblad operators acting on the system Lq\ we need to create a set of 
size \Sp\ Lindblad operators - one to operate on each subspace independently: 

8(L<§) = {L<§®m\,jeS P }. (A.4) 

At first glance it might appear that we have increased a system of size m — \S\ to 
size mn, where n = \Sp\. Whilst this is true, the situation is not as bad as it seems 
at first, because the coherences between the different subsystems are unimportant - 
instead of a density matrix of size (nm) 2 we can use a system of size nm 2 , an increase 
linear in the number of system states. In practice we can often do better than this by 
eliminating unnecessary states from some of the subsystems. 

As a simple concrete example, consider a resonantly driven two-level system 
Sq — {\g) , |e)}, with Hamiltonian 

H Q = Sl(\g)(e\ + \e)(g\). (A.5) 
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Suppose that it is also possible for the system to spontaneously emit from state \e) into 
the environment - a transition described by the jump operator 

J Q = \g}(e\. (A.6) 

We are interested in knowing about the time distribution of the first time a photon 
is emitted. We add the two process-states and 1, indicating whether the event has 
occurred or not. Our new system is described by: 

S = {\g0) , |e0) , \gl) , |el)} (A.7) 

H = fi(|<70)(e0| + |e0)(<?0| + |^l)(el| + |el)( 5 l|) (A.8) 

Ji = |<?l)(eO| (A.9) 

J 2 = |<7l)(el|. (A.10) 

The population in the 1 subspace at time t will give us the probability a photon has 
been emitted by this time. 

If we only cared about this distribution, we could reduce the size of the system by 
eliminating the state | el) , removing the last two terms from the Hamiltonian as well as 
J 2 . 
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